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(i) 


ABSTRACT 


This  thesis  deals  with  stability  criteria  and  Lyapunov 
functions  for  certain  second  order  differential  equations. 

The  function  x"  +  x  =  f(t,x,e)  is  investigated  for 
stability  of  its  periodic  motions  by  transforming  the  problem  to  an 
equivalent  stability  problem  of  equilibrium  points. 

In  Chapter  I,  some  properties  such  as  singular  points  and 
criteria  for  their  classification  are  reviewed.  Definitions  are  given, 
Lyapunov  functions  are  defined  and  Lyapunov  theorems  are  stated. 

Chapter  II  is  devoted  to  the  transformed  Duffing  equation. 

The  stability  of  this  equation  is  examined. 

In  Chapter  III,  the  transformed  Van  der  Pol  equation  is 
investigated.  Here  reference  is  made  to  the  existence  of  a  limit  cycle. 
Both  the  phase  plane  and  the  Lienard  plane  are  discussed. 

In  Chapter  IV,  the  transformed  simple  pendulum  equation 
is  discussed  first  without  damping  and  then  with  damping. 

Finally  some  of  the  effects  of  perturbation  on  the  stability 
of  a  system  are  discussed  in  an  Appendix  at  the  end  of  the  thesis. 
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CHAPTER  I 


BASIC  CONCEPTS 


In  this  chapter  we  review  topics  pertaining  to  the  development 
of  the  other  chapters.  Reference  is  made  to  specific  texts  where 
further  details  can  be  found. 


1.  Study  of  Singular  Points. 

Consider  a  second  order  system  of  the  form 


dx 

dt 


P(x,y) 


Q(x,y). 


(1.1) 


Such  a  system  in  which  the  independent  variable  does  not  appear 
explicitly  is  called  autonomous. 

We  shall  assume  that  P(x,y)  and  Q(x,y)  are  defined  in 
some  domain  D  of  the  x  -  y  plane  and  satisfy  a  Lipschitz 
condition  in  both  x  and  y  in  some  neighbourhood  of  each  point 


Hurewicz  [6]  has  shown  that  if  tQ  is  any  real  number  and 


of  D. 


■ 


2. 


(xQ,y0)  any  point  of  D,  there  exists  a  unique  solution  of  (1.1) 

x  =  x ( t )  y  =  y ( t)  (1.2) 

satisfying  x(t  )  -  x  y(t  )  =  y 

o  o  J  o  ^o 

Furthermore  it  is  known  that  either 

(i)  The  solutions  (1.2)  are  defined  for  all  real  values  of  t  or 

(ii)  If  the  solution  (1.2)  are  not  defined  for  t  >  t^  (say),  then 
either 

(a)  As  t  ->  t^-0,  the  point  (x(t),y(t))  approaches  the  boundary 
of  D  or  if  D  is  unbounded,  possibly 

(b)  As  t  -»■  t^-0  either  x(t)  or  y(t)  or  both  become  unbounded. 

In  any  subdomain  D*  of  D  in  which  P(x,y)  does  not 
vanish  we  may  write  (1.1)  in  the  form 


=  Q(x>y)  , 

dx  P(x,y) 


(1.3) 


This  is  a  restricted  direction  field  and  we  know  that  through  any 
point  of  D’,  there  exists  a  unique  integral  curve  of  (1.3),  for 
if  two  distinct  integral  curves  were  to  have  a  point  in  common,  they 
would  have  to  be  tangent  at  this  point,  a  possibility  ruled  out  by 
the  fact  that  P  and  Q  are  Lipschitzian.  Let  us  consider  a 
solution  (1.2)  of  (1.1)  with  the  restriction  that  x(t)  and 
y(t)  are  both  constants.  Then  (1.2)  defines  a  curve  with  continuous 


3. 


turning  tangent.  This  curve  is  called  a  characteristic  curve  or  a 
trajectory.  Note  there  is  a  difference  between  solutions  and 
trajectories  of  (Id)  :  a  trajectory  is  a  curve  in  D  that  is 
represented  parametrically  by  more  than  one  solution.  Thus  x(t), 
y(t) ,  x(ti-c),  y(t+c)  represent  distinct  solutions  but  they  represent 
the  same  curve  parametrically. 

We  may  think  of  (Id)  as  a  "flow"  in  the  phase  plane 
defined  by  a  velocity  vector  field.  Consider  the  vector  field 
defined  by  V(x,y)  -  (P(x,y) ,Q(x,y) )  that  is  at  (x,y),  the  vector 
V(x,y)  has  a  horizontal  component  P(x,y)  and  a  vertical  component 
Q(x,y).  Equation  (Id)  defines  the  motion  of  a  particle  (x,y)  in 
the  plane  by  the  condition  that  its  velocity  at  each  point  be  equal 
to  the  vector  V(x,y). 

The  above  material  is  fully  discussed  in  Hurewicz  [6] . 

Definition  1.4,  A  singular  or  a  critical  point  ^x0’yc^  a 
stationary  point  of  the  flow  PCx^y^)  -  Q(xo>yo)  =  0  and  the 
integral  curves  passing  through  it  consists  of  the  point  itself. 

Definition  1.5  A  critical  point  (x_,y_)  of  (1.1)  is  called 

-  Q  O 

an  isolated  critical  point  if  there  exists  a  neighbourhood  of  (x^y^) 
which  contains  no  other  critical  points. 

We  now  introduce  the  notion  of  stability  of  a  critical 
point  on  equivalently  stability  of  the  solution  x(t)  -  xq,  y(t)  -  VQ 
— oo  <  t  <  05  of  (Id)  . 


' 


1 


-  ri-’  -  \1. 


4. 


Definition  1.6 ,  Let  ( x  ,y )  be  an  isolated  critical  point  of 
(1.1),  then  (xo>yo)  is  said  to  be  stable  if  given  any  e  >  0,  there 
exists  6  >  0  such  that 

(i)  every  trajectory  of  (1.1)  in  the  6-neighbourhood  of  (x^y^) 
for  some  t  =  t^  is  defined  for  t^  t  <  °°  and 

(ii)  if  a  trajectory  satisfies  (i)  it  remains  in  the  e-neighbourhood 

of  (x  ,y  )  for  t  >  t,  . 
o  o  1 

If  in  addition  every  trajectory  C  :  x  =  x(t),  y  35  y(t) 

satisfying  (i)  and  (ii)  also  satisfies  (iii)  lim  x(t)  =  x  and 

t  -*  » 

lim  y(t)  =  y  ,  then  (x  ,y  )  is  said  to  be  asymtotically  stable. 

t  ->  oo 

Finally,  an  isolated  critical  point  which  is  not  stable  is  said  to  be 
unstable i 


The  definition  of  stability  roughly  states  that  (xo,yQ)  is 
stable  if  once  a  trajectory  enters  a  small  disc  containing  (xQ,yo) 
it  remains  within  a  larger  disc  for  all  future  time.  The  above 
definition  is  sometimes  called  stability  to  the  right.  A  similar 
definition  is  given  for  stability  to  the  left  as  t  . 

2 .  Singularities  of  a  linear  system. 

In  this  section  we  consider  the  linear  system 

dx  i  i 

^ “  ax  +  by 
=  cx  +  dy  . 


(1.7) 


. 


5. 


where  a,b,c,d  are  real  constants.  Therefore  we  may  let  D  be  the 
entire  xy-plane  and  so  all  solutions  are  uniquely  defined  on 
-oo  <  t  <  „  Hence  we  can  discuss  the  behaviour  of  trajectories  in 

the  phase  plane  of  (1.7).  We  discuss  a  system  of  the  form  (1.7) 
because  a  complete  description  of  the  phase  plane  can  be  given. 
Secondly  many  systems  can  be  expressed  in  the  form 


x*  -  ax  +  by  +  e^(x,y) 
y'  -  cx  +  dy  + 

(the  prime  (’)  denotes  differentiation  with  respect  to  t) .  If 
and  are  sufficiently  small  in  a  neighbourhood  of  a  critical  point, 

we  hope  chat  the  behaviour  of  the  trajectories  is  locally  like  that 
of  (1.7).  The  system  (1.7)  has  a  singularity  at  (0,0).  Assume 
it  has  no  other  singularity  that  is 

|  a  b  j 

i  |  4  o 
i  c  d  | 


The  characteristic  equation  of  (1.7)  is 


2 

“  A 


(a+d)X  +  (ad-bc)  =  0 


c 


d  -  X 


' 


where  A  -  ~2~  —  ~2  ^ ^a~^2  +  •  Several  cases  arise  depending  on 

2 

the  value  of  the  discriminant,  A  ~  (a-d)  +  4bc.  These  cases  will  be 

discussed  in  a  manner  similar  to  that  of  Ritzer  and  Rose  [9]. 

Case  I .  A  =  0.  This  could  occur  if  a  =  d  ^  0  and  b  =  c  =  0. 
Thus  x'  =  ax  j  y'  s  ay 


x 


at  at 

c^e  y  = 


Therefore  c^x  -  c^y„ 

When  all  the  solutions  go  through  the  singularity  as  in  this  case  we 
have  a  node.  If  a  <  0  as  t  -+■  »,  x  and  y  -*■  0,  the  node  is  stable. 
If  a  >  0,  it  is  unstable o 

Case  II. 


A  >  0  b  ~  c  «  0  ,  a  ^  d 

We  get 

x '  ~  ax  y'  =  dy 

at  dt 

x  -  c^e  y  -  c^e  . 


Eliminating  t,  we  get 


' 


7. 


a  d 

V  "  C3X  ^ 

If  a  and  d  have  the  same  sign  we  have  a  node.  If  a  and  d  are 
of  different  signs  so  that  when  x  -»■  0,  y  -*  «>,  when  x  «>,  y  0, 

we  have  a  saddle  point.  Whether  we  have  a  node  or  saddle  point,  there 

are  two  special  cases  which  are  obtained  when  c^  -  0  or  when  c^  =  0. 

The  curves  are  then  y  ®  0  and  x  =  0.  In  the  node,  the  curve  x  =  0 

is  the  only  one  which  enters  the  origin  vertically.  In  the  saddle  the 
curves  x  =  0  and  y  =  0  are  the  only  ones  which  enter  the  singularity. 

Case  III.  A  <  0.  This  occurs  if  d  =  a  and  c  -  -b  .  The 
differential  equation  becomes 

x '  -  ax  +  by 

y ’  =  -bx  +  ay  * 

These  equations  have  a  solution 

at  .  ,  at  ,  . 

x  -  e  sin  bt,  y  -  e  cos  bt. 


Hence 


If  a  <  0,  the  distance  /x2  +  y^  from  the  origin  to  a 


8. 


point  moving  on  the  solution  curve  will  decrease  steadily  as  the  point 
moves  round  the  origin.  In  such  a  case,  the  singularity  is  called  a 
stable  spiral  curve. 

If  a  >  0,  the  curves  spiral  away  from  the  origin  and  we  have 
an  unstable  spiral. 

In  the  above  example,  if  a  ~  0,  the  solution  would  be 
x  =  c^  sin  bt  y  =  c^  cos  bt.  The  integral  curves  will  be  circles 
on  more  generally  ellipses.  In  this  case  the  origin  is  a  centre. 


Table  1  ;  Classification  of  Singularities 


^2  = 
dx 


cx-tdy 

ax-rby 


or 


(  x '  -  ax  +  by 
]  y»  s  cx  +  dy 


f  '  2 

Fundamental  quantities  jf  A  -  (a-d)  +  4bc 

<  p  s  a  +  d 


ad  -  be  ^  0 


Case 


Singularity 


Stability 


(1)  A  =  0 

(2)  A  >  0 

(3)  A  <  0 


Node 


j Node  if  q  >  0 
Y  Saddle  if  q  <  0 

j  Centre  if  p  =  0 
^Spiral  if  p  ^  0 


Stable  if  p  <  0 
Unstable  if  p  >  0 

Stable  if  p  <  0 
Unstable  if  p  >  0 


j  Stable  if  p  <  0 
\  Unstable  if  p  >  0  . 
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3 .  Stability  for  Non-Autonomous  Systems. 

Consider 


x'  =  f(t,x)  (1.8) 

where  x  =  x(t)  -  (x. (t)  ,  .  .  .  ,x  (t))  is  an  unknown  n-dimensional 

i  n 

vector  function  and  assume  that  f(t,x)  =  (f . (t ,x) , . . . ,f  (t,x)) 

1  n 

is  defined  and  continuous  in 

D  -  1  (t  ,x)  |  r1  <  t  <  «,  ||  x  1 1  <  aj  . 

-  ■ 

By  j  ( x 1 1  we  mean  jjx||  ~  (  l  j  x  j  )  .  A  solution  (not  necessarily 

i=l 

unique)  of  (1,8)  satisfying  x(t^)  “  XQ  will  be  denoted  by 

x  ( t )  *  x  (  t  j  t  ,  x  )  . 

o  o 

Definition  1.9  Let  x(t)  =  x(t,t  ,x  )  be  a  solution  of  (1.8) 
-  o  o 

satisfying 

(i)  x(t)  is  defined  for  t  <_  t  <'  00  and 

(ii)  the  point  (t,x(t))  belongs  to  D  for  t  >  t  . 

Then  x(t)  is  said  to  be  stable  if: 

(a)  There  exists  y  >  0  such  that  every  solution  x(t,'to,x^) 

satisfies  (i)  and  (ii)  ,  whenever  jjx^-x^lj  <  y  an^ 

(b)  given  e  >  0,  there  exists  6>0,  0  <  6  _<  y  such  that  |  j  x  Q— x  ^  |  j  <  6  implies. 

||  x  ( t }  t  j  x  )— x(t,t  jX,)  1 1  <  £  for  t  <  t  <  ». 
o  o  ox  o 


A  solution  that  is  not  stable  is  said  to  be  unstable. 


10. 


Definition  1.10.  The  solution  x(t)  =  x(t,t  ,x  )  of  (1.8) 
-  o  o 

Is  asymtotically  stable  if  it  is  stable  and  in  addition  there  exists 

P  >  0,  0  <  p  <_  y  such  that  |jxQ-XjJj  <  p  implies  lim  j  j  x  ( t ;  t  ,xq) 

t  -*■  » 

-  x(t;to,x1)  ||  =  0. 

Geometrically,  the  definitions  say  that  x(t)  is  stable  if 
any  other  solution  whose  initial  data  is  sufficiently  close  to  that 
of  x(t)  remains  in  a  "tube"  enclosing  x(t).  If  the  diameter  of 
the  tube  approaches  zero  as  t  -»•  00 ,  then  x(t)  is  asymtotically 
stable . 

4=  Lyapunov’s  Direct  or  Second  Method. 

We  now  discuss  an  important  method  of  studying  the  equation 
(1.8)  .  The  method  is  known  as  Lyapunov’s  direct  or  second  method 
and  involves  constructing  a  particular  type  of  function  V(t,x)  from 
which  the  stability  or  instability  of  the  solution  in  question  can  be 
determined. 

Assume  f(t,x)  =  (f . (t ,x)  ,  ,  .  .  ,f  (t,x))  satisfies  the  following 

1  n 

conditions 

(i)  f(t,x)  is  defined  and  continuous  in  D  -  {(t,x)  |  ||xjj  <  a, 

r^  <  t  <  «•'} 

(ii)  a  condition  assuming  uniqueness  of  solutions  x(t:to>x^)  is 
satisfied  at  every  point  (tQ,xo)  in  D  and 

(iii)  f(t,0)  *  0  for  all  t  and  hence  x(t;t  ,0)  ~  0  is  a  solution 

o 

of  (1.8)  for  t  >  r. 

o 


. 
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Definition  1.11,  The  class  K  consists  of  all  continuous  real 
valued  strictly  increasing  functions  <?>(r),  0  <_  r  <_  a  ,  which  vanish 
at  r  -  0  . 


Let  0  <  b  £  a  and  t  >  r^,  suppose  V(t,x)  is  a  real 
valued  function,  continuous  together  with  its  first  partial 
derivatives  in  the  set 

B  =  {  (t  ,x)  |  t  £  t  <  CD,  |  j  x ||  £  b} 

Furthermore  assume  V(t,0)  =  0  for  t  >  t  . 

—  o 

Definition  1.12  The  function  V(t,x)  is  positive  (or  negative) 
definite  if  there  exists  a  function  <p  in  the  class  K  such  that 

V(t,x)  _>  <p  (|  |  x|  j )  or  V(t,x)  <_  — <p  (|  |x|j  )  for  all  (t,x)  in  B. 

Definition  1.13,  The  function  V(t,x)  is  said  to  be 
descrescent  or  to  admit  an  infinitesimal  upper  bound  if  there  exist 
h  >  0  and  a  function  tj;  in  the  class  K  such  that,  jv(t,x)j  £  ^(|(x||) 
for  |jx||  -  h  and  t  £  t  . 

By  V*  we  denote  the  function 


V’ 


V'(t,x) 


n 

V 

L 

i=l 


d  v 

dx  , 

1 


f^^CtjX) 


We  will  now  state  some  theorems  on  stability,  the  proofs  of 


' 

t 
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these  may  be  found  in  Sanchez  [10]. 

Theorem  1 . 14,  If  a  function  V(t,x)  exists  that  is  positive 
definite  and  whose  derivative  V'  is  non  positive,  then  the  solution 
x(t)  h  0  of  (1.8)  is  stable. 

Theorem  1.15  If  V(t,x)  exists  that  is  positive  definite, 
descrescent  and  whose  derivative  is  negative  definite,  then  the 
solution  x(t)  =  0  of  (1.8)  is  asymtotically  stable. 

Definition  1016  A  function  V(t,x)  satisfying  the  hypotheses 
of  Theorem  1.14,  is  called  a  Lyapunov  function  of  equation  (1.8). 

In  the  case  of  the  autonomous  system,  we  omit  the  dependence 
of  V  on  t  and  delete  the  family  K  which  serves  to  ensure 
uniformity  with  respect  to  t.  Thus  V  =  V(x)  is  a  positive 
(negative)  definite  function  if 

(i)  V  is  continuous  together  with  its  partial  derivatives  in  some 
neighbourhood  of  the  origin  and 

(ii)  V(x)  >_  0  or  (<_  0)  with  equality  when  x  =  0, 

The  two  theorems  in  the  autonomous  case  are  much  easier  to 
prove.  Now  we  will  be  dealing  with  the  case  n  =  2.  The  curve 
V(x,y)  =  c  constitutes  a  family  of  concentric  loops  encircling  the 
origin.  The  ovals  grow  in  size  as  c  increases  from  zero  and  each 
oval  encloses  others  of  smaller  c. 


The  hypotheses  of  Theorem  1«14  now  imply  that  given  a  critical 


' 
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point  at  the  origin  and  a  suitable  Lyapunov  function  in  an  open 
region  R  about  it s  all  solutions  that  start  within  the  oval 
will  remain  in  it. 

Theorem  1.15  now  implies  that  if  5,  V(x,y)  =  c,  does  not 
cross  the  boundary  of  R,  all  solution  trajectories  that  start  at  a 
point  P^Cx^jy^)  will  tend  to  the  origin  as  t  ^ . 

We  have  spoken  only  of  stability.  A  natural  question  to  ask 
is,  "When  do  we  have  instability?" 

In  the  autonomous  case  we  just  need  V  to  be  positive  definite 
V*  to  be  positive  definite  and  in  every  neighbourhood  N  of  the 
origin  V(xq)  >  0  for  x^  e  N.  This  gives  instability. 

5.  We  shall  discuss  the  method  we  shall  use  in  examining  the  stability 
of  the  Duffing,  Van  der  Pol  and  simple  pendulum  equations  (the  equations 
with  which  we  will  be  dealing) . 

Consider  again  equation  (1.8)  .  If  the  unperturbed  motion 
is  defined  as  the  solution  p(t,t  ,xq)  of  (1.8),  then  its  stability 
can  always  be  reduced  to  the  stability  of  the  equilibrium.  For  this 
purpose,  we  set  x  =  y  +  p(t,to,xQ)  in  (1.8)  and  obtain,  the 
differential  equation 

y ’  -  f (y+p(t,to,xQ) , t)  -  f (p(t,tQ,xo) ,t)  (1.17) 

The  right  side  of  this  equation  is  a  function  of  y  and  t  which 


t 
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vanishes  for  y  =  0.  The  stability  behaviour  of  the  solution 
p(t,to,xo)  is  equivalent  to  the  stability  behaviour  of  the  equilibrium 
of  (1.17).  Equation  (1.17)  is  called  the  differential  equation  of 
the  perturbed  motion. 

The  goal  of  this  transformation  is  a  simplification  of  the 
stability  problem;  however  the  new  differential  equation  is  often  more 
complicated  than  the  original  one.  So  the  differential  equation  of  the 
perturbed  motion  for  the  solution  of  an  autonomous  differential 
equation  is  almost  always  non-autonomous . 

In  what  follows,  the  new  equation  will  be  considered  on  its 
own  merit.  Certain  conditions  will  be  imposed  on  p(t)  to  get 
suitable  Lyapunov  functions. 


CHAPTER  II 


THE  TRANSFORMED  DUFFING  EQUATION 


1. 


In  engineering  equations  of  the  type 


mx"  +  ij>(x')  +  f(x)  -  f(t) 


(2.1) 


are  studied.  The  term  mx"  is  frequently  referred  to  as  the  inertia 
force,  -(^(x1)  as  the  damping  force,  -f(x)  as  the  restoring  force 
or  spring  force  and  f(t)  as  the  external  force.  The  terms  in 


equation  (2.1)  are  interpreted  in  amny  ways  in  the  case  of 


electrical  and  combined  electrical-mechanical  systems  of  various  types. 


The  equation 


x"  r  x  +  &x^  =  f (t) 


(2.2) 


is  an  equation  of  this  type  and  is  called  the  Duffing  equation 
(without  damping) .  It  is  one  of  our  purposes  in  this  chapter  to 
study  the  stability  of  motions  of  this  equation.  This  is  done  by 
means  of  the  transformation  discussed  in  the  last  chapter. 

Let  p(t)  be  a  solution  of  (2.2)  .  We  use  the 
transformation  y  ~  x  -  p(t),  giving  x  =  y  +  p(t) . 


16. 


Hence  xf  =  y'  +  p'(t)  ,  x"  -  y"  +  p"(t)  .  Substituting 
into  equation  (2.2),  we  have 

y"  +  p"(t)  +  y  +  p(t)  +  £(y3+3y2p(t)+3yp2(t)+p3(t) )  =  f(t). 


Therefore 


y"  +  y  +  3  (y3i-3y2p(t)+3yp2(t)  ) 

=  -  p"(t)  -  6p3(t)  -  p(t)  +  f(t)  =  0.  (2.3) 

Now  y  h  0  is  a  solution  of  equation  (2.3)  corresponding  to  x  =  p(t) 
and  hence  the  problem  of  stability  is  reduced  to  studying  stability 
of  equilibrium  points. 

It  is  therefore  natural  to  investigate  the  stability  of  the 
origin  of  the  transformed  equation,  including  those  conditions  chat 
should  be  put  on  p(t)  to  obtain  results  concerning  the  nature  of 
the  stability. 

Actually,  we  will  be  trying  to  find  suitable  Lyapunov 
functions . 

Without  loss  of  generality,  let  us  consider  the  equation 
x"  +  x  +  3(x3+3x2p(t)+3xp2(t))  =  0 


(2.4) 


U)q 
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instead  of  equation  (2.3). 

Consider  first  the  simpliest  form  of  this  equation,  the 
case  when  p(t)  =  0.  Then  equation  (2.4)  becomes 

x"  +  x  +  gx3  =  0  .  (2.5) 

In  the  phase  plane,  this  becomes 

x '  =  y 

y’  -  -x  -  gx3  (2.6) 

Now,  a  Lyapunov  function  can  be  thought  of  as  an  "energy"  function. 

A  suitable  Lyapunov  function  for  equation  (2.6)  is 

V(x ,y)  =  -  y  +  x  +  -^  x  ■*  h.  (2.7) 

There  are  two  cases  (non  linear)  to  consider,  g  *  0  and  g  <  0, 

referred  to  as  the  hard  and  soft  spring  case  respectivley .  For  the 
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case  g  =  0,  V(x,y)  =  —  x  +  —  y  is  a  Lyapunov  function.  In 
this  case  we  have  stability  in  the  whole  plane  —  this  shall  be 
referred  to  this  thesis  as  global  stability. 

Hard  Spring  g  >  0, 


The  origin  is  the  only  critical  point  of  the  system.  V(x,y) 


18. 


given  by  (2,7)  is  the  Lyapunov  function  we  use, 


We  again  have  in  this  case  global  stability.  The  trajectories  in 
the  phase  plane  are  ovals  symmetric  about  the  origin.  The  linear 
approximation 


shows  that  the  origin  is  a  centre  and  the  ovals  are  traversed  counter 
clockwise.  This  is  discussed  in  Brand  [2] . 

Soft  Spring  g  <  0; 

Equation  (2.6)  has  three  critical  points  (0,0)  and 
(0,+  x  )  where  x^  =  /-1/6  . 

The  origin  is  still  a  centre  but  we  do  not  have  stability 
in  the  whole  plane  (global  stability)  for  the  points  Pq  =  (0,x^) 
and  =  (0,-x^)  are  now  saddle  points.  If  we  shift  the  origin  to 
Pq  by  putting  x  =  X  +  x^  and  y  =  Y  we  get 


X'  -  Y 


Y'  =  -  (X+xx)  -  3(X+X;l) 


3 


. 


' 
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whose  linear  approximation  is 


X '  -  Y 

Y'  -  -  (l+33x^)X  =  2X • 


From  the  table  in  Chapter  1,  we  see  that  this  is  a  saddle  point. 


2 .  Curves  of  constant  energy 


We  will  discuss  curves  of  constant  energy  in  the  (x,y)  plane. 


The  closed  energy  curves  are  integral  curves  corresponding  to  periodic 
motion  x(t).  If  x(t)  is  periodic,  the  corresponding  x  -  y  curve 
is  closed.  Conversely  if  any  x  -  y  curve  is  closed,  it  follows  that 
the  displacement  and  velocity  at  any  time  t  are  reached  again  after 
a  certain  time  T,  that  is,  x(t+T)  =  x(t) .  Hence  the  motion  is 


periodic . 


In  a  neighbourhood  of  the  origin  the  curves  given  by  equation 
(2.7)  are  closed  curves  which  have  the  appearance  of  ovals  since 


3  —  can  be  neglected  in  comparison  with 


1  2 

—  x  for  small  x.  See 


Stoker  [11]. 


In  the  case  of  the  soft  spring  3  <  0 


v  2  2  . 2  x  2 

V(x,y)  =  y  +  x  -  k  -y  -  h  ,  k  =  -  3* 


2 


(2.8) 


to  J  arirdei  ■  j  sn  i.  * n*  "  v  .  .j  [  -*  ■  3  ms  as.!  rc  ) 
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2 

The  right-hand  side  of  equation  (2.8)  is  a  quadratic  in  x  with 

2  2 

discriminent  1  -  2k  h;  it  is  therefore  positive  if  1  -  2k  h  is 

negative  and  in  this  case  h  must  of  necessity  be  positive.  The 

transition  from  curves  which  cross  the  x-axis  to  curves  which  do 

2  1 
not,  occurs  at  1  -  2k  h  =  0  that  is  at  h  =  — j  . 

2k 


On  this  curve  the  velocity  y^  corresponding  to  x  =  0 
1  1 

the  value  y  =  -  .  So  for  y  >  -  we  have  open  curves,  y  < 

°  k/2  °  k/2 

we  have  closed  curves  encircling  the  origin.  The  transition  curve 

corresponding  to  y  =  has  two  double  points  on  the  x-axis  at 

°  k/2 

x  =  +  /-1/3  =  +  /l/k^  =  1/k  .  This  curve  separates  the  x  -  y 

plane  into  regions  in  which  three  different  types  of  curves  occur. 
See  Figure  2.1  below. 
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Figure  2.1. 
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Figure  2.1  shows  two  sets  of  curves  which  represent  non¬ 
periodic  motion  and  one  set  of  curves  which  represent  periodic 
motion. 

We  can  show  that  the  separatrices  for  the  soft  spring  are 
parabolas.  The  separatrices  are  the  "dotted"  curves  in  the  figure. 


A(y  -  ^  )  =  x^ 

k/2 


A  = 


o  u  /  /2  1  \  2 

So  we  have  ( - r-  y  -  — r)  -  x  or 

R  kZ 


k  ,  !  2.. 


-k  1  2 

The  other  parabola  is  y  =  —  (— j  -  x  )  .  The  region  of  stability  is 

/ 2  k 

the  region  bounded  by  these  two  curves . 


3 .  Case  p(t)  constant, 


Equation  (2.4)  becomes 
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(2.9) 


Lyapunov  function,  V'(x,y)  =  0.  Thus  for  the  case  3  >  0,  we 
again  have  stability,  the  origin  is  a  centre.  When  3  <  0,  we  claim 
that  we  have  two  saddle  points  as  in  the  case  p(t)  =0.  We  will 
also  examine  the  region  of  stability  if  one  exists.  One  critical 
point  is  the  origin,  we  can  find  the  other  two  critical  points  by 
setting  y'  =  0. 


Then 


3x  +  3p3x  +  3p  3x  +  x  =  0 


3  ,  „  „  2  ,  „  2 


x(3x^+3p3x+3p^3+l)  =  0 


x 


2 


2 


Hence  the  other  two  critical  points  are 


, 
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P  = 


(-  f p  i 


-3  2  1 


p  -  j  ,  0)  and 


<-  f  p  +, 


'-3  2 

T  p 


?  •  0) 


Here  we  must  stipulate  that 


3  2  1. 

IP  - ? ' 0 


i.e. 


3  2  1 

Ip  "e 


or 


2 

3P 

4 


1 

6 


3p‘ 


<  0  <  0  . 


Let  us  examine  the  system  (2.9)  by  putting  0 
We  find  that 


m  ,4  .. 

—  (3  >  m  >  0) 


x  =  y 


y*  a  _  x  +  JL  (x^+3px2+3p^x) 
P 


3  2 

a  _  x  +  ™L_  +  3mx_  +  3mx 

P2  P 


3  2 

,  „  , s  ,  mx  ,  3mx 

-  (3m-l)x  +  — 7T  +  - 

P 

P 


. 


( 
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Examining 


0  1 
(3m-l)  0 


the  linear  part  of  the  system,  we  see  that 

A  -  4 (3m-l) 
q  =  -  (3m-l) 


1-4  1 

when  A  >  0  i.e.  m>-^-  ( — j  <  3  < - j) ,  we  have  a  saddle  point 

3P^  3pZ  x 

at  the  origin.  When  m  <  1/3  (0  >  3  >  -  — j)  we  have  a  centre  at 

x  3pZ 

the  origin.  For  m  >  — 


P  “  (-  j  P 
o  z 


f  3 
("  2  p 


('3p ,  0). 


Thus  (~  |  P  -  /-  -  "2  +  2~ 


4  *  m 


3  2 

"  4  P 

2 

+  £- 
m 

/  3  2 

„  2 

"  4  P 

+  3p 

=  (~3p 

+  6  , 
m  ’ 

)  .  Similarly 

m 


pi  ■  <-  i  p  + 


1 p2  +  >  0) 


=  (-«  ,0). 

m 


“ic'i 


£*  k  -  V-  if'1  *  o'5 
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When 


m 


< 


1 

3 


0) 


=  (-3p  -  6m,  0) > 


2 

P 


~ T 

+  £- 
m 


0)  «  (6  ,0). 

m 


We  can  now  examine  the  singularities  at  each  of  these  points.  The 

point  (-6  ,0)  terms  out  to  be  a  centre  and  the  other  three  points 
m 

(-3p+6  ,0),  (-3p-6  ,0)  and  (6  ,0)  are  all  saddle  points, 
mm  m 


In  conclusion  we  can  say  that  the  regions  of  stability  is 
found  wherever  we  have  a  centre.  The  two  figures  given  below  should 
be  useful. 

Y  Y 

A  A 


X 


m 


X 


Figure  2.2, 


Figure  2.3 


0  >  3  >  - 

3p 
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In  the  light  of  the  preceeding  discussion  we  get  the  following  theorem. 

Theorem  2.12.  V(x,y)  =  —  y  +  -j  x  +  3  ("^  x  +px  +  y  P  x  )  is  a 

suitable  Lyapunov  function  for  the  system  (2.9)  with  3  >  0  .  For 

-4  1 

each  value  of  p  and  for  — y  <  3  < - y  ,  we  have  a  saddle  point  at 

3p  3pZ 

the  origin  of  the  corresponding  system.  The  other  two  singularities  of 

1 

the  system  are  a  centre  and  a  saddle.  For  0  >  3  >  -  — y  ,  we  have  a 

3p 

centre  at  the  origin,  the  other  two  singularities  are  saddles. 


4 .  The  case  p(t)  non  constant. 

x'  =  y 

y'  =  -(l+33p2(t))x  -  3(x3+3p(t)x2)  . 

2 

Let  A(t)  =  1  +  33p  (t),  we  will  first  examine  the  system 


x 


t 


“  y 


y’  -  -  (A ( t ) ) x 


(2.14) 


with  A(t)  >  0  and  then  with  A(t)  <  0  .  A(t)  >  0  implies 

2  2  i 

(l+33p  (t))  ■>  6  >  0  .  Equation  (2.14)  can  now  be  written  as 


x"  +  (l+33p2 (t) )x  =  0  . 


Multiplying  by  x' 


dx 

dt 


we  get 


j  fj  X  11 J  4  9W  «  (3)  '*  ^  '■  +  ; 
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x"  +  Cl+3Bp2(t))x  i|  =  0 


Integrating  we  get 


(^2~)  2  +  /  (l+3Bp2 (t) )xx’ dt  =  . 

But  (l+3Bp2(t))  5  2  and  x*  =  y  .  Therefore 


O^)  '  +  /  62xx'dt  <_  (-^)2  +  f  (i+3Bp2(t))xx'dt  = 

2  2 

+  5  —  C1  (2.15) 

It  is  clear  that  for  x  and  y  satisfying  (2.15)  we  have  ellipses 
surrounding  the  origin.  In  this  case  the  origin  is  a  centre  and  we 
have  stability. 

When  A(t)  <  0,  ( l-»-3 3p 2  (t) )  -  62  <  0  .  We  then  get 


^r- +  /  (l+SBp^  (t)  )xx  ’  dt  =  C2 


V  -  6x2 
—  4 


(2.16) 


Equation  (2.16)  represents  hyperbolas  and  we  have  a  saddle  point  at 
the  origin  —  an  unstable  situation.  We  state  the  following  simple 


. 


. 
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theorem. 

Theorem  2.17-  Given  system  (2.14),  we  have  stability  at  the  origin 
if  A(t)  >  0  and  instability  if  A(t)  <  0. 

The  system  (2.13)  as  a  whole  will  now  be  studied,  we 
will  try  to  find  a  suitable  Lyapunov  function  that  will  give  us  some 
information  concerning  the  nature  of  the  stability.  Consider 


Tr/  ,  12,12,3*  2 ,  ,  2,84 

V(x,y,t)  =  2  x  +  2  y  *  2  +  T  X  ‘  (2.18) 

(We  first  look  at  the  case  8  <  0) . 

Now  (i)  V(x,y,t)  is  defined  for  all  t  >_  0 

(ii)  V(0,0,t)  =0  for  t  >_  0 

(iii)  V (x  ,y ,  t)  >_  W(x,y)  . 

W(x,y)  will  now  be  specified.  Let  0  <  p(t)  <  M.  Then 

tT/  1  2  ,  1  2  ,  3  oyj2  2  ,  8  4  ,0  10. 

V(x,y,c)  >  y  x  +7y  +  2  x  +‘4X*  (2.19) 

W(x,y)  is  to  be  positive  definite.  Thus  it  should  be  that 

(1+38M2)  >  0 
38  M2  >  -1 
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The  right  hand  side  of  (2.19)  is  of  the  form 


Ax2  +  ^  y2  +  3  where  A  =  (1+3  3M) 


and  since  this  is  to  be  our  W(x,y),  it  is  necessary  that 


.  2  .  .  x 

Ax  +  3  ~t~  y  '0 
4 


x2(A  +  S  0 


x_  _  A 

4  <  3 


2 

x  < 


4A 

3 


'-4A 


<  x  < 


-4A 


So  for  the  case  3  <  0,  the  region  in  which  V(x,y,t)  is  defined 


must  be  such  that 


-H 


(1+33M2) 


fl 


(1+33M  ) 


For  stability  we  want  V’(x,y,t)  <_  0 


V'(x,y,t)  -  xx'  +  yy*  +  3Sp^(t)xx’  +  33p ' (t)p (t)x2  +  3x2x' 


ft 
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=  xy  -  xy  -  33p2(t)xy  -  Sx3y  -  33p(t)x2y 
+  33p2(t)xy  +  33p(t)p' (t)x2  +  3x3y 
=  33x2(p' (t)p(t)-p(t)y) 

We  want  (p’(t)p(t)  -  p(t)y)  _>  0 

p’(t)  >  y  . 

We  will  show  that  under  the  assumption  0  <  p(t)  <  M  the  origin  is 
unstable.  For  stability  at  the  origin  we  should  look  at  a  neighbour¬ 
hood  of  the  origin  say  |xj<6,  | y |  <6  .  Thus  we  have 

P*  (t)  _>  6 

This  implies  that  p(t)  becomes  unbounded.  So  for  negative  3  we 
are  unable  to  find  appropriate  assumptions  on  p(t)  under  which  the 
origin  is  stable. 

For  3  >  0  ,  the  positive  definite  function  is 
12  12 

W(x,y)  =  —  x  +  —  y  .  Note  in  this  case  the  positive  definite  function 
was  easy  to  find  in  the  sense  that  no  restrictions  had  to  be  placed 
on  x  . 


We  will  consider  p(t)  _>  6  >  0  .  For  stability 
p*  (t)p(t)  -  p (t)y  <_  0 


i.e . 


p'(t)  1  y 
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Since  we  are  concerned  with  stability  at  the  origin,  for  y  _<  6  , 
p'(t)  <_  6  and  thus  contradicts  the  fact  that  p(t)  >  0  .  If 
p(t)  <_  -6  <  0  ,  for  stability 

p'  (t)  _>  y  . 

Thus  p'  (t)  _>  6  and  contradicts  p(t)  <  0  .  We  can  also  see  that 
for  positive  £  ,  but  are  unable  to  find  appropriate  assumptions  on 
p(t)  under  which  the  origin  is  stable. 


' 


CHAPTER  III 


THE  TRANSFORMED  VAN  DER  POL  EQUATION 


1 .  Historical  Note. 

The  Van  der  Pol  equation 

x"  +  y(x^-l)x'  +  x  =  0 

arose  in  the  study  of  oscillatory  vacuum  tube  circuits.  Van  der  Pol 
showed  that  certain  properties  of  vacuum  tubes  could  be  predicted 
only  if  the  non  linear  terms  of  the  equation  were  retained.  The 
linear  theory  predicted  that  no  periodic  oscillations  in  the  current 
would  occur  under  certain  conditions  when  experimental  evidence  showed 
they  did  occur.  By  retaining  non-linear  terms  in  the  equation,  Van 
der  Pol  was  able  to  give  a  complete  explanation  of  the  phenomenon. 

2.  We  will  be  considering  the  equation 


x"  4-  yCx^-^x1  -I-  x  =  f(t). 


(3.1) 


Let  x  =  p(t)  be  a  solution.  Using  the  tranformation  y  =  x  -  p(t) 
or  x  =  y  +  p(t)  we  get  x’  =  y'  +  p'(t)  and  x"  -  y"  +  p"(t). 
Substituting  into  equation  (3.1)  we  have 
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y"  +  p"(t)  +  p (y2+2p(t)y+p2(t)-l) (y'+p' (t))  +  y  +  p(t)  =  f(t)  (3.2) 
Now  since  p(t)  is  a  solution 

p"(t)  +  p (p2(t)-l)p’ (t)  +  p(t)-  f(t)  =  0  . 

Hence  we  get  the  resulting  equation 

y"  +  u (y2+2p(t)y+p2(t)-l)y'  +  (up ' (t)y+2pp (t)p ' (t)+l)y  =  0  . 

This  is  the  equation  we  will  study.  Without  loss  of  generality,  let 
us  consider 

x"  +  u (x2+2p (t)x+p2 (t)-l)x'  +  (up* ( t)x+2pp (t )p ’ (t)+l)x  =  0.  (3.3) 

First  we  consider  the  case  when  p(t)  3  0,  then  p(t)  ~  p,  a  constant 
and  lastly  p(t)  as  a  non-constant  function  of  t.  Before  doing  this, 
some  material  pertaining  to  the  understanding  of  this  chapter  will  be 
introduced . 

3.  We  consider  the  concept  of  a  limiting  set.  This  was  first  introduced 
by  G.D.  Birkoff. 

Let 


x*  =  X(x)  ,  X(0)  =  0  . 


(3.4) 
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Intuitively,  if  x(t)  is  a  solution  of  equation  (3.4),  its  positive 
limiting  set  I  is  whatever  the  curve  x(t)  tends  to  with  infinite 
time.  Thus  if  x(t)  spirals  around  6  ,  a  limit  cycle,  5  is  its 
positive  limiting  set.  (In  the  plane  a  limit  cycle  is  an  isolated 
closed  path  corresponding  to  a  periodic  solution) .  If  it  tends  to  a 
point  A,  the  point  is  its  positive  limiting  set. 

+ 

More  accurately  p  is  in  F  if  there  is  an  increasing 

sequence  of  times  t  -*  °°  with  n  such  that  x(t  )  -*•  p  as  n  +  <*. 

n  n 

If  x(t)  is  bounded,  then  x(t)  approaches  its  positive  limiting  set 
r+  as  t  -*  a*  ,  that  is,  given  any  e  >  0,  if  N(e)  is  the 

-f- 

e-neighbourhood  of  F  ,  then  there  is  a  time  T  such  that  for  t  >  T, 
x(t)  lies  in  N(e).  This  definition  is  due  to  LaSalle  and  Lefschetz 
[7].  We  now  give  another  definition. 

Definition  3.5  An  invariant  set  G  is  characterized  by  the 
property  that  if  a  point  xq  is  in  G  then  its  whole  path  (forward 
and  backward)  lies  in  G. 

From  the  definition  of  limiting  set  there  follow  also  the 
property:  if  x(t)  is  bounded  for  t  >_  0  and  if  a  set  M  contain 

F+,  then  x(t)  tends  to  M  as  t  «>. 

4.  The  Van  der  Pol  equation  is  a  special  case  of  the  equation 


x"  +  f(x)x’  +  g(x)  »  0 


, 
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which  is  known  as  the  Lienard  equation.  This  equation  has  been  shown 
to  have  a  unique  stable  limit  cycle  under  the  following  conditions; 

(a)  f(x)  even,  g(x)  odd,  both  continuous  for  all  x  and 

f (0)  <  0. 

(b)  xg(x)  >  0  for  x  f  0. 

(c)  For  every  interval  |x|  <  k,  there  is  an  h  such  that 

|g(x  )-g(x  ) |  <  hjx^-x^l,  a  Lipschitz  condition 
x 

(d)  F(x)  =  J  f(x)dx  as  x  00 . 

o 

(e)  F(x)  has  a  single  zero  at  x  =  a  and  is  monotone  increasing 
for  x  >_  0.  For  a  proof  of  this  theorem  one  can  see  Cesari  [4]  . 

LaSalle  and  Lefshetz  [7]  have  proven  two  theorems  which  deal 
with  the  extent  of  asymtotic  suability  of  the  Lienard  equation.  We 
state  these  theorems. 

Theorem  3 . 6  Lee  V(x)  be  a  scalar  function  with  continuous 
first  partial  derivatives.  Let  ft  designate  the  region  where 

X/ 

V(x)  <  l  .  Assume  that  ft  is  bounded  and  that  within  ft 

X,  X/ 

(a)  V(x)  >  0  for  x  0 

(b)  V'  (x)  <_  0  0 

Let  R  be  the  set  of  all  point  within  ft^  where  V’(x)  =  0.  Let  M 
be  the  largest  invariant  set  in  R.  Then  every  solution  x(t)  in 
ft.  tend  to  M  as  t  ■> 

Xj 


ll  ;  >i3oJmvea  Jo  x  &  U 
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Theorem  3.7.  If  condition  (b)  of  the  above  theorem  is  replaced 
by  V’(x)  <  0  for  all  x  #  0  in  A  then  the  origin  is  asymtotically 

Xj 

stable  and  above  all  every  solution  in  A,  tends  to  the  origin  as 

t  -*■  00. 

By  considering  the  function 

V(x,y)  -  \  y2  +  G(x) 


as  a  Lyapunov  function  for  the  Lienard  equation  and  putting  conditions 
on  G(x) ,x ,g (x)F (x) ,  they  showed  that  the  region  A^  defined  by  V  <  i 
was  a  measure  of  asymtotic  stability. 

The  Lienard  equation 

x"  +  f(x)x*  +  g(x)  =  0 


is  equivalent  to 


x"  =  y  -  F (x)  ,  y*  =  -g(x) 

for  differentiating  we  get 

x"  =  y*  -  F'(x)x*  =  g (x)  -  f(x)x’. 

When  y  =  x’  +  F(x) ,  the  x  -  y  plane  is  called  the  Lienard  plane. 


. 

' 
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Except  when  f(x)  =  0,  it  differs  from  the  phase  plane  in  which  y  =  x' . 


5.  We  return  to  equation  (3.3)  and  the  case  p(t)  =0.  We  now  get 

x"  +  y(x2-l)x’  +  x  =  0. 


This  is  the  Lienard  equation  in  which 


f(x)  =  y(x2-l)  ,  F(x)  =  y(j  x3-x) 


g(x)  -  x  . 


The  stability  of  the  solutions  depends  on  the  sign  of  y  . 


Case  1.  y  >  0  .  The  phase  plane  equations 


x  ’  =  y 

2 

y'  -  -x  -  y (x  -l)y 


(3.8) 


have  a  critical  point  (0,0).  When  x  <  1 


V(x,y) 


2 


-  x 


2 

+  y 


is  a  Lyapunov  function,  for 


' 
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V'(x,y)  =  2(xx'+yy’)  =  -2y(x2-l)y2  <  0  . 


2  2 

The  interior  of  the  circle  x  +  y  =1  is  a  region  of  asymtotic 
stability.  Any  trajectory  that  starts  at  a  point  within  the  circle  tend 
to  the  origin.  The  linear  part  of  (3.8)  has  the  matrix 


I  0  1  ’ 

l-l  M  ■ 

From  the  table  in  Chapter  I, 

A  2 
A  —  y  -4 

p  =  p 

q  =  1  . 

So  for 


y  <  -  2  the 
y  =  -  2  the 
-  2  <  y  <  0  the 


origin  is  a  stable 
origin  is  a  stable 
origin  is  a  stable 


node 

node  (inflected) 
spiral . 


Case  II.  y  >  0 

Now  f(0)  =  -y  <  0 .  The  conditions  for  the  existence  of  a 
stable  limit  cycle  are  satisfied.  The  origin  is  now  unstable.  The 
-trajectories  for  any  postive  y  are  same  for  -y 


reversed  in  sense. 
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All  trajectories  inside  the  limit  cycle  tend  to  merge  with  it  as  t  -*■  ». 

2  2 

The  limit  cycle  encloses  the  circle  x  +  y  =1  in  the  phase  plane. 

We  will  now  consider  the  situation  in  the  Lienard  plane, 


t  ,13., 

x  =  y  -  yU  x  -x)  y  = 


-x 


12  2 

If  y  <  0,  V(x,y)  =  (x  -t-y  )  is  a  Lyapunov  function  for 


V ' (xy)  =  yx2(l-  -j  x3)  <  0 


2  2  2 
when  x  <  3  .  All  paths  within  the  circle  x  +  y  =3  tend  towards 

the  origin.  For  corresponding  positive  jy|  ,  all  such  paths  are  reversed 

in  sense,  cross  the  circle  and  eventually  merge  with  the  limit  cycle. 


For  all  positive  values  of  y  >  the  limit  cycle  of  the 

Van  der  Pol  equation  in  the  Lienard  plane  enclose  the  circle 

2  .  2  , 
x  +  y  =  3 . 


Brand  [2]  discusses  the  Lienard  equation. 


6 .  The  case  p(t)  =  p  (a  constant). 

Equation  (3.3)  now  has  the  form 


x"  +  y ( ( x+p ) 2 — 1 ) x  *  +  x  =  0 


x  slosio  9ri3  alrfJiw  ari  ; bo  XIA 


' 
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or 


x'  =  y 

y'  =  -y ( (x+p) 2-l)y  -  x  .  (3.9) 

We  want  to  study  stability  of  this  system. 

Let  V(x ,y)  =  ■—  x2  +  y2  .  Then 

V '  (x  ,y)  =  xx  '  +  yy  ’ 

=  xy  -  y  ( (x+p)  2-l)y2  -  xy 
=  -  y ( (x+p) 2  -  1)  . 

2 

For  y  <  0  ,  V*  <_  0  if  (x+p)  <_  1  .  This  implies  -1  _<  x  +  p  _<  1 
thus  -p-l_^x_<l-p  .  Now  we  are  concerned  with  stability  at  the 
origin,  therefore  we  would  like  x  to  be  in  an  interval  containing 
the  origin.  For  this  to  happen,  p  has  to  be  suitably  chosen. 

Choose  p  such  that  -1  <_  p  _<  1  ,  that  is  j  p  |  _<_  1 
Therefore  the  region,  in  which  we  will  be  considering  the  Lyapunov 
function,  is  bounded  on  the  right  by  the  line  x  =  1  -  p  and  on  the 
left  by  x  =  -1  -  p  . 

2 

We  get  instability  if  (x+p)  >  1 

i.e.  x  +  p  >  1  ,  x  >  1  -  p 


' 
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or  x  +  p  <  -1  ,  x  <  -1  -  p  . 

Y 


Figure  3.1 

The  shaded  area  is  the  region  where  the  Lyapunov  function  is  considered. 

When  y  >  0  ,  V1  (x,y)  <_  0  if  (x+p)^  _>  1  ,  that  is 

(a)  x  _>  1  -  p  or 

(b)  x  _<  ~1  -  p 

Now  we  are  concerned  with  stability  at  the  origin,  so  in  case  (a)  we 
choose  p  such  that  p  _>  1  and  in  case  (b)  p  <_  1  .  We  thus  get 
as  show  on  the  following  page. 
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Y 


Figure  3.2 

2 

We  have  instability  if  (x+p)  <  1 


Y 


i .  e . 


-1  <  x  +  p  <  1 


-l-p<x<l-p 


7 .  The  case  p(t)  a  non  constant  function  of  t  . 

Equation  (3.3)  can  now  be  written  as 


x'  =  y 

y'  =  -y( (x+p(t) )2-l)y  -  (pp * (t)x+2yp(t)p ' (t)+l)x  . 


Consider  the  function, 


V (x,y , t)  =  \  x2  +  \  y2  +  "3  P'  (t)x3  +  yp (t) p '  (t)x2 


-1  -  p 

- 


(3.10) 


(3.11) 


* 
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We  shall  show  that  this  is  a  Lyapunov  function  for  the  system  (3.10) 
by  imposing  certain  conditions  on  p(t)  .  We  consider  here  only 
the  case  p  >0.  We  claim  that 

(a)  V(x,y,t)  is  defined  for  all  t  >_  0  in  some  region. 

(b)  V(0,0,t)  =  0  for  t  >_  0. 

(c)  V(x,y,t)  >_W(x,y)  (a  positive  definite  function). 

Condition  (b)  is  obviously  true  so  we  have  now  to  find  the  region 
and  also  W(x ,y) . 

Now  V(x,y,t)  =  y2  +  x2(-|  +  j  p'  (t)x  +  yp(t)p’  (t))  . 

Let  us  define  U(t,x,y)  such  that 


U(t,x,y)  =  j  +  y  p'(t)x  +  yp(t)p'(t)  . 


If  we  choose,  p(t)  >_  0,  0  <_  p’(t)  <  M,  and  restrict  x  to  be 

2  2 

in  the  region  -a  <  x,  a  >0,  then 


U(t,x,p)  >_  j  p’  ( t ) x 


1  y  w  2 

>  2  “  3  Ma 


(Note  since  we  are  concerned  with  stability  at  the  origin  we  want  x 


9f(i  :'0 


. 
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to  be  in  some  interval  containing  the  origin) . 


2  1 

Furthermore  if  we  choose  0  <  6  <  —  ,  we  get  that 

1  U  „  2  r2  2  3/2  -  362 

Ma  >6  if  a  <  - — -  . 

2  3  —  —  yM 


With  these  conditions  imposed. 


,  N  12  2  2 
W(x,y)  =  j  y  +  6  x 


We  finally  have  to  show  that  V(x,y,t)  <_  0,  then  expression  (3.11) 
will  be  a  Lyapunov  function  and  the  origin  will  be  stable. 


V'(x,y,t)  =  xx’  +  yy’  +  ^  p"  (t)x3  +  yp’(t)x2x*  +  yx2  (p  ( t)  p'' ( t)  +  (p  '  ( t)  )  2) 

+  2yp(t)p' (t)xx’ 

2  2  2  2  2  2  2 
=  xy  -  ux  y  -  2yp (t)xy  -  yp  (t)y  +  yy  -  yp'(t)x  y 

-  2yp(t)p* (t)xy  -  xy  +  y  p"(t)x3  +  yp'(t)x2y 
+  yx2(p(t)p"(t)+(p’ (t))2)  +  2yp(t)p' (t)xy 

=  yx2(p(t)p"(t)+(p’(t))2)  -  y(p2(t)-l)y2 

-  y (x2y2+2p(t)xy2-p" (t)  y-  ). 


(3.13) 


.  .  ■  ' 
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2  2  2  x 

We  will  consider  -y (x  y  +2p(t)xy  - —  p"(t))  as  higher  order  terms. 

Equation  (3.13)  is  of  the  form 


V’(x,y,t)  =  A(t,y)x2  -  B(t,y)y2  -  y (x2y2+2p(t)xy2-  ~  p'r(t)) 


where  A(t,y)  =  y (p(t)p" (t)+(p 1 (t) ) 2) ,  B(t,y)  =y(p2(t)-l).  Since 
y  >  0,  for  stability,  it  is  necessary  to  have  A(t,y)  <_  0,  B(t,y)  >_  0 

Thus  p(t)p"(t)  +  (p'(t))2  £  0 


(Note  that  we  needed  p(t)  to  be  twice  continuously  differentiable) 

2 

i.e.  p(t)  e  C  . 

2 

B(t,y)  >_  0  implies  p  (t)  1  therefore  p(t)  _<  -1  or 

p(t)  >2  1  •  Since  we  are  considering  p(t)  to  be  positive,  we  take 
only  the  latter  condition  and  replace  the  original  condition  on  p(t), 
0  <_  p(t)  ,  by  1  <_  p(t)  . 

With  these  conditions  imposed  on  p(t),  we  see  that 
expression  (3.11)  is  a  Lyapunov  function  for  the  system  (3.10), 
for  the  higher  order  terms  can  be  neglected  as  will  be  shown  presently 

Since  we  are  considering  stability  at  the  origin,  we  should 
take  into  account  values  of  x  and  y  in  a  small  region  around  the 
origin.  Our  main  purpose  is  to  find  out  what  effect,  if  any,  do  these 
higher  order  terms  have  on  the  nature  of  the  stability.  We  claim  that 


1  •  '• 

1  't  O  i 

■  : vx(:  )a£-fV‘x)M  •  v(u-  )•  -  -  (3 
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they  can  be  neglected  under  certain  conditions. 

2  2 

The  term  -yx  y  obviously  causes  no  trouble  for  it  is 

2  x^ 

always  negative  for  positive  y  .  We  will  look  at  -2p(t)xy  +  —  p"(t). 
Consider  |x|  <  e^,  |  y  [  <  £2  ei  ’  £2  >  ^  and  very  small,  let 
e  -  max  (e^^e),  we  now  have  that 


2  £ 
-2x(p(t)y2+  p"(t))  £  -2e1(p(t)£2+  —  p"(t)) 


<  ke' 


(k  a  constant) 


3  2  2 

Now  ke  is  small  compared  with  A(t,y)x  -  B(t,y)y  which  is  of  the 

2 

form  k^e  provided  e  is  sufficiently  small. 

In  the  light  of  the  previous  discussion,  we  can  now  state 
the  following  theorem. 


Theorem  3.14t 


(i) 

p(t)  e  C2,  1 

(ii) 

0  <  p*  (t)  <  M 

(iii) 

.2  1 

0  <  6  <  j  > 

(iv) 

2  2 

x  >  -a  ,  a  £ 

yM 


The  function  V(x,y,t)  defined  by  expression  (3.11)  is  a  suitable 
Lyapunov  function  for  the  system  (3.10). 


A  typical  p(t)  satisfying  the  conditions  of  the  theorem 


'?■  ■' '  ’3  *‘rfT 

' 

ovsri  won  ®w  ,(3  .  a)  x  m 

3  ^ 
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is  p(t)  =  3  -  e 

We  should  note  that  we  can  find  many  functions  p(t) 
which  satisfy  the  conditions  of  Theorem  3.14.  There  also  exists 
various  Lyapunov  functions  for  the  system  (3.10)  and  these  are 
obtained  by  imposing  the  necessary  conditions  on  p(t). 

Although  the  case  y  <  0  will  not  be  dealt  with  here, 
it  is  possible  to  find  a  Lyapunov  function  for  th&  case  and  formulate 
a  theorem  similar  to  Theorem  3.10. 

Also,  the  system  (3.10)  ,  really  points  out  the  fact 
that  the  "new"  differential  equation  formed  by  the  transformation 
is  more  complicated  than  the  original.  This  can  be  seen  by 

comparing  system  (3.10)  and  system  (3.8). 

It  should  be  noted  too,  that  in  finding  the  Lyapunov 
function  for  this  nonautonomous  case,  the  positive  definite 
function  W(x,y)  played  a  role,  whether  indirectly,  in 
determining  the  region  of  stability. 

Finally,  it  should  be  realised  that  there  may  be 
values  of  y  ,  positive  or  negative,  in  this  non-autonomous 
case  and  a  function  V  which  could  present  an  unstable 
situation  at  the  origin.  It  is  possible  that  we  could  have 
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stability  only  in  a  region  bounded  away  from  the  origin. 


We  will  now  discuss  the  simple  pendulum  equation. 


CHAPTER  IV 


THE  TRANSFORMED  PENDULUM  EQUATION 


1.  A  simple  pendulum  consists  of  a  mass  m  supported  by  a  weightless 
rod  of  length  £  pivoted  at  the  top.  The  position  of  the  mass  is  given 
by  the  angle  x  that  the  rod  makes  with  the  vertical.  The  equation  of 
motion  is  given  by 


x"  +  a  sin  x  =  0  (4.1) 

were  a  =  y  and  g  is  the  acceleration  due  to  gravity. 

In  this  chapter,  we  will  consider  equation  (4.1)  and  two 
others  then  make  a  transformation  and  study  the  resulting  equations. 
The  second  equation  is 

x"  +  ex'  +  a  sin  x  =  0  .  (4.2) 

This  is  the  damped  simple  pendulum  equation  and  the  last  equation  is 
the  forced  damped  simple  pendulum  equation 

x"  +  ex'  +  a  sin  x  -  ef(t)  .  (4.3) 

In  these  three  equations  a  >  0,  c  >  0  and  in  equation  (4.3), 
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no.  tfiDps  ujlJbnsq  alqnrxs  b9qm£.b  b£>  vie  :  \ 


'  xo  +.  ' 
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f ( t+T )  =  f(t). 

Before  going  into  the  discussion  of  these  three  equations, 
we  will  discuss  nonlinear  conservative  systems.  This  material  can 
be  found  in  Minorsky  [8],  also  in  Butenin  [3]. 

2 .  Fundamental  properties  of  non-linear  conservative  systems. 

The  simpliest  differential  equation  for  such  a  system  is 

x  "  +  f  (x )  =  0  (4.4) 


which  is  equivalent  to 


x '  =  y  (4.5) 

y'  «  -f(x) 

and  the  differential  equation  of  the  integral  curves  is 

(4.6) 

dx  y 

which  shows  that  the  integral  curves  have  a  horizontal  tangent  at  points 

x  which  are  roots  of  f(x)  =  0,  provided  y  ^  0  at  these  points.  As 
i 

to  singular  points,  they  require  simultaneously  y  =  0,  f(x)  =  0.  The 
energy  integral  is 


nlflsJuk  «i  o ,,  .b  «  >  o  1  J  K, 


i  -  >  . 
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~  y2  +  n(x)  =  h 


(4.7) 


kinetic  energy  and  n(x)  as  the  potential  energy,  h  the  total  energy 

expresses  the  fact  that  the  system  is  conservative.  For  a  given  h, 

equation  (4.7)  represents  a  trajectory  in  the  phase  plane  which  exists 

as  long  as  h  -  IT  (x)  >0.  A  position  of  equilibrium  corresponds  to  the 

case  y  =  0,  n '  (x)  =  0.  The  latter  condition  implies  the  differential 

equation  has  an  extremum  at  the  equilibrium  point.  Now  equation  (4.7) 

1  2 

is  equivalent  to  y  y  =  h  -  II  (x)  .  Let  us  examine  the  auxiliary  xz 
plane  on  which  we  shall  construct  z  =  n (x)  ,  (see  Figure  4.1)  and  draw 
the  line  z  =  h.  Rea<l  motion  takes  place  only  for  those  x  located 
to  the  left  of  the  noint  of  intersection  of  z  =  II  (x)  and  the  line 
z  =  h  .  Now  v/v/2  =  +  -  n(x)  ,  it  follows  on  the  phase  plane 

(x ,y ^  =  y//2)  the  integral  curves  will  have  the  form  shown  in 


Figure  4.2, 


z 


n(x) 


7\ 


h  -  n(x) 


h 


1 


>  X 


Figure  4.1. 


, 
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A 


\ 


, _ \ 

+/h-n  (x)  \ 


\ 


-Vh-Ib 


/ 


>  X 


Figure  4.2 


The  Duffing  equation  discussed  in  Chapter  II  is  an  example 
of  a  conservative  system. 


As  discussed  in  Minorsky  [8],  the  stability  of  equilibrium 


is  associated  with  the  extremum  values  of  n(x)  ,  the  potential  energy. 
When  n (x)  is  a  minimum,  the  curves  in  the  neighbourhood  of  the  origin 
are  ellipse  so  the  origin  is  a  centre  and  when  II  (x)  is  a  maximum, 
the  curves  in  the  neighbourhood  of  the  origin  are  hyperbolas,  the  origin 
is  a  saddle  point. 

3.  Returning  to  equation  (4.1)  and  letting  x  =  y  +  p(t)  we  get 


y"  +  p"(t)  +  a sin(y+p (t) )  =  0 


y"  -  a  sin  p(t)  +  a  sin(y+p(t))  =  0 


. 

•  JUS  ' 8V8  7.  .  1  .  .  '■ 

, 


" 
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y"  +  a(sin(y+p(t) )-sin  p(t))  =  0 


y"  +  a(2  cos  ^ )  sin  ^  )  =  0 


Without  loss  of  generality,  we  consider 


x"  +  2a  cos  (X+2P  ^)sin  ^  =  0 


(4.8) 


In  this  case 


ix  = 

dx 


-2a  cos  ^ )  sin  y 


We  have  a  conservative  system  here.  When  p(t)  =  0,  we  get 


dy  _  a  sin  x 
dx  y 


ydy  =  -  a  sin  xdx 

2 

•^r-  =  h  -  a  /  sin  x  dx 


JI(x)  =  a  /  sin  x  dx  =  -  a  cos  x  . 

2 

Therefore  y  /2  =  h  +  a  cos  x  ,  n (x)  has  isolated  minima  for 
x  =  0 ,  +  2ir ,  +  4tt  , .  . .  .  The  minima  correspond  to  stable  equilibrium 
states  in  the  phase  plane.  They  are  centres.  When  x  =  +tt,+  3tt,...  , 

we  have  isolated  maxima  (saddle  points). 
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For  h  <  -a,  real  motion  of  the  pendulum  does  not  exist. 

For  -  a  <  h  <  a  ,  the  trajectories  are  closed  and  surround 
singular  points  which  are  centres.  Closed  trajectories  correspond  to 
periodic  solutions  of  the  pendulum. 

For  h  =  a  ,  the  trajectories  are  separatrices  which  intersect 
at  saddle  points.  For  h  >  a,  the  trajectories  are  divergent.  But 
the  x-coordinate  is  periodic,  the  motion  of  the  representative  point 
along  the  trajectories  correspond  to  a  clockwise  or  counter-clockwise 
motion  of  the  pendulum  about  the  suspension  point.  When  p(t)  =  nn, 

X  X 

for  n  even,  cos  (y  +  p(t))  =  cos  y  and  the  singularities  are  the  same 

X  X 

as  before.  For  n  odd  cos  (y  +  p(t))  =  -  cos  y  •  In  this  case  we 
have  maxima  at  0 , 2  tt,4 tt,  .  . .  ,  and  hence  saddle  point  where  we  had  centre 
and  vice  versa. 

4.  The  case  when  p(t)  =  p,  a  constant,, 

Equation  (4.8)  can  be  written  as 

XXX 

x"  +  2a(cos  —  cos  p-sin  —  sin  p)sin  y  =  0 

2  x 

x"  +  acos  p  sin  x  -  2a  sin  y  sin  p  =  0 

x"  +  a  cos  p  sin  x  -  a(l-cos  x)sin  p  =  0 


x"  +  a (cos  p  sin  x+sin  p  cos  x)  -  a  sin  p  =  0. 


jJ  or.  qI-J  j!  :  j.  >  jjU 

. 
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Or  x '  =  y 


y*  =  -a(cos  p  sin  x  +  sin  p  cos  x  -  sin  p) 


dy  __  -g(cos  p  sin  x  +  sin  p  cos  x  -  sin  p) 
dx  y 


i.e.  ydy  =  (-a  cos  p  sin  x  -  a  sin  p  cos  x  +  a  sin  p)dx 

2 

y 

=  a  cos  p  cos  x  -  a  sin  p  sin  x  +  ax  sin  p  . 


In  this  case 


II(x)  =  -  J  -  a  cos  p  sin  x  -  a  sin  p  cos  x  +  a  sin  p  dx 


=  -a  cos  p  cos  x  4a  sin  p  sin  x  -ax  sin  p 

-  -acos(p4x)  -ax  sin  p  . 

A  position  of  equilibrium  corresponds  to  the  case  when 
H’(x)  =  0,  that  is  when 

sin  (x+p)  -  sin  p  =  0 


and  so 


p  +  x  =  p  4  2Ttn 


- 
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x  =  2nir  , 


Or  when 

p  +  x  =  7T  -  (p-27in) 
x  =  (2n+l)TT  -  2p 


i .  e . 


x  =  rnr  -  2p  (n  odd)  . 

Now  n"(x)  =acos  (x+p)  and  when  n"(x)  >  0,  we  have  a  minimum  (a 
centre);  for  n"(x)  <  0,  we  have  a  maximum  (a  saddle  point). 

If 

x  =  2nn,  n"(x)  =acos  (p+x)  =acos  p  . 


If 


x  =  utt  -  2p,  n"(x)  =acos  (p+x)  =acos  (p+nTT-2p) 

=acos  (nTT-p) 

=  (-1)  acos  p 


=  -acos  p  . 


So  if  cos  p  >  0,  we  have  centres  at  x  =  2un  for  all  n.  At  this 


time  we  have  saddle  points  at  mr  -  2p  for  n  odd.  The  converse 
happens  when  cos  p  <  0. 
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5 .  The  case  p(t)  a  non  constant  function  of  t „ 

Recall 

x"  +  2a (cos  (.*+2P(t))  sin  f)  =  0. 
Dr 


(4.8) 


x '  =  y 

y'  =  -2a(cos  (^2P(t))  slnf).  (4.9) 

We  consider  the  function, 

12  2 

V(x,y,  t)  =  I  y  +a  CGS  P(c)  (4.10) 

This  will  be  shown  to  be  a  Lyapunov  function  for  the  system  (4.9) 
with  suitable  restrictions  on  p(t) . 

Firstly, 

(i)  V(x,y,t)  is  defined  for  t  >_  0.. 

(ii)  V(0,0,t)  =  0  for  t  0  . 

That  these  two  conditions  hold  is  obvious.  We  further  claim  that 

12  x2 

V(x,y,t)  dominates  W(x,y)  =  —  y  +  ak  —  where  k  will  be  specified 
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2 

yy*  +  a  cos  p(t)xx'  -  a  -y-  p’(t)  sin  p(t) 

-i  2 

-LX  X 

-2ay  cos  —  (x+2p(t))sin  y  +  a  cos  p(t)xy  -a  y-  p'(t)sin  p(t) 

2 

-ay  [  sin(x-fp  (t)  )-sin  p(t)]  +  axy  cos  p(t)  -ay-  p'(t)sin  p(t) 

-ay[sin  x  cos  p(t)+cos  x  sin  p(t)-sin  p(t)]  +  axy  cos  p(t) 

2 

-  ay-  p’  (t)  sin  p(t) 


-ay  cos  p(t)[sin  x  -  x]  -  ay  sin  p(t)[cos  x-1] 
2 

-ot  -y-  p*  (t)  sin  p(t) 


x  5  x  4 

-ay  cos  p(t) (-  yy  +  0(x  ))  -  ay  sin  p ( t) [ -  y —  +  o(x  ;] 


x 


-  a  —  p' (t)  sin  p(t) 


y  x^y  cos  p(t)  -  axZ  sin  p(t)  (— y— —  yO  -  ay(cos  p(t)o(x‘)) 


2  . 


T  (t) 


4. 

sin  p(t)o(x  )) 


2 

2,-axy  cos  p(t)  ,  ap'(t)  sin  p(t)N  ,  x  .  /i_N 

-  x  ( - 1 -  y - +  a  -y-  y  sin  p(t) 


5  4 

-  ay(cos  p(t)o(x  )+sin  p(t)o(x  ))  . 


' 
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2 

This  is  of  the  form  -x  (-B (t)xy+A( t) )  +  higher  order  terms. 

Here  A(t)  =  fn  ?<*>  ,  B(t)  =  g-££g-P^  and  the  hlgher  order 

terms  are 


2  5  4 

ax  y  sin  p(t)  -  ay(cos  p(t)o(x  )  +  sin  p(t)o(x  )) 


For  stability,  it  is  necessary  that 


A(t)  -  B ( t)xy  >  0 


i .  e . 


xy 


A(t) 

B(t) 


We  also  impose  the  condition  that 


A(t)B(t)  >  0 


for  if  not  we  have  xy  <  n,  n  negative  and  since  we  are  considering 
stability  at  the  origin,  we  would  not  want  this  to  occur. 

ap’(t)sin  p(t)  a  cos  p(t)  s 

Thus  — ^ — - - — ^ -  •  - -j—*- -  >  0  .  We  choose  p(t) 

such  that 

0<6<sinp(t)<_l-6*  0<6<1,0<6'  <1 
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Then  0  <  arc  sin  6  £  p(t)  <_  arc  sin  (1  -  6') 

0  <  cos  arc  sin  (1  -  6’)  £  cos  p(t)  £  cos  arc  sin  6  <  1  . 

(Note  we  can  now  see  that  the  k  we  need  in  our  W(x,y)  is 

12  2 

k  =  cos  arc  sin  (1  -  6 ’ ) .  Thus  V(x,y , t)  £  —  y  +  —  cos  arc  sin  (1  -  6  '  )x 
We  also  choose  p’(t)  such  that 

q  >  p ’ (t)  >  0 


(q  will  be  chosen  suitably  small)  and 


Then 


lim  p '  ( t)  =  0  . 

t  -*•  OO 


xy 


ctp’  (t)  sin  p(t) 

2 

■g  cos  p(t) 


=  3p ’ (t)  tan  p(t) . 


Now 


cos  arc  sin  6  — 


tan  p(t)  < 


1-6’ 


—  cos  arc  sin(l  -  6’) 


So 


0  <  3p ’ (t) tan  p(t)  £ 


3g(l  -  6’) 

✓26'  -  6’2 
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Therefore 

we  get  that 


xy  < 


3q (1  -  6’) 

/2 6’  -6’2 


If  we  choose  q  =  /2s*  - 


xy  <3(1  -  6')  =  £ 


£  >  0 


With  this  restriction  on  xy  we  have  stability  at  the  origin. 


Y 


> 


X 


Figure  4.3, 

(Shaded  portion  shows  stable  region.) 

It  remains  to  discuss  the  higher  order  terms.  We  will 

show  that  these  are  easily  dispersed  of  when  we  are  close  to  the  origin. 

5  4 

It  is  obvious  that  in  the  higher  order  terms,  -ay(cos  p(t)o(x  )  +  sin  p(t)o(x  )) 

ii  ii  2 

can  be  neglected  for  |x|  <  ,  |y|  <  •  The  term  ax  y  sin  p(t)  is 

3  2 

of  the  form  Ce  which  is  small  when  compared  with  -x  (A(t)-B(t)xy) 

2 

which  is  of  the  form  C^e  provided  e  is  sufficiently  small, e  =  max (£^,£2)* 
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In  the  light  of  the  preceeding  discussion,  we  state  the  following 

theorem. 

Theorem  4.11  Given  (i)  p(t)  e  C', 

(ii)  0  <  6  £  sin  p(t)  £  1  -  6'  0  <  6  <1,  0  <  6’  <  1 

(iii)  0  <  p'(t)  <  q,  lim  p’(t)  =  0  . 

t  ->  00 

The  function  V(x,y,t)  defined  by  (4.10)  is  a  suitable  Lyapunov  function 
for  the  system  (4.9). 

6.  Performing  the  usual  transformation  on  equation  (4.2)  and 

changing  y's  for  x's  we  get 

x"  +  ex'  +  2  (cos  (x+— )  sin  y)  =  0  .  (4.12) 

We  are  considering  here  a  dissipative  system  for  we  have  large  deviations 
of  the  pendulum  subject  to  friction. 

Cape  p(t)  s  0.  Equation  (4.12)  can  now  be  written  as 

x '  »  y 

y’  =  -cy  -  a  sin  x  •  (4.13) 

In  the  phase  plane  ■  -c  -  a  -Syn— -  .  The  singularities  are  x  =  kv , 
y  >=  0  (k  is  any  integer  positive  or  negative)  and  correspond  to 
stable  or  unstable  equilibria  according  as  k  is  even  or  odd.  Taking 
the  first  term  in  the  expansion  of  sin  x  ,  the  system  (4.13) 


becomes 


, 
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x 


? 


■  y 


y '  =  -ax  -  cy  . 


2 

Here  A=c  -4a,  p  =  -c,  q=a.  Referring  to  Table  1,  Chapter  I, 

2 

we  see  that  if  c  >  4a,  we  have  a  stable  node  at  the  origin.  If 
2 

c  <  4a,  we  have  a  stable  spiral. 

The  singularities  which  were  centres  in  the  case  where  there 
was  no  damping  become  spiral  or  nodal  points.  We  saw  in  the  earlier 
case  that  at  x  =  tt  ,  the  singularity  was  a  saddle  point.  In  the  case 
with  damping,  we  again  get  a  saddle  point  for  the  singularity  is  not 
affected  by  the  damping  term  ex'  .  The  singularities  at  x  =  kfT  are 
stable  if  k  is  even,  unstable  if  k  is  odd.  The  tangents  to  the  paths 
are  vertical  on  the  x-axis  and  horizontal  on  the  curve 


y 


a  sin  x 
c 


On  the  basis  of  these  facts,  one  can  construct  a  qualitative  graph  of 
the  arrangement  of  the  trajectories.  From  the  figures  below,  it  follows 


Figure  4.4  Stable  spiral 


Figure  4.5  Stable  node 
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It  is  interesting  to  note  that  part  of  the  phase  plane  between 
-tt  and  it  can  be  thought  of  as  laid  out  on  the  surface  of  a  cylinder. 

We  will  briefly  discuss  the  cylindrical  phase  space.  This 
is  fully  dealt  with  in  Andronow  and  Chaikin  [1], 

The  description  of  the  behaviour  of  a  dynamical  system  by 
means  of  a  phase  space  requires  a  one-one  correspondence  between  the 
states  of  the  system  and  the  points  of  the  space.  Up  to  now,  we  have 
been  investigating  physical  systems  whose  phase  space  can  be  represented 
by  a  plane.  In  general,  however,  this  plane  will  not  be  satisfactory, 
as  is  seen  by  considering  the  system  of  the  ordinary  pendulum  (the 
system  with  which  we  are  dealing) .  Here  the  state  of  the  system  is 
determined  by  the  angle  of  deviation  from  the  equilibrium  position  and 
the  velocity  of  the  pendulum.  If  we  restrict  ourselves  to  motion  which 
does  not  exceed  one  revolution,  no  difficulties  will  arise  if  the  phase 
space  is  taken  as  a  plane.  But  when  the  angle  of  deviation  is  increased 
by  a  multiple  of  2tt  ,  we  obtain  a  state  identical  with  the  original 
state.  In  the  phase  "plane",  we  would  have  an  unlimited  number  of  points 
corresponding  to  the  same  physical  state  of  the  system.  Since  the 
necessary  one-one  correspondence  does  not  hold,  a  plane  is  not  in 
general  a  suitable  phase  space  for  a  pendulum.  In  particular,  for 
systems  whose  position  is  completely  defined  by  their  deviation  angle, 
the  condition  of  one-one  correspondence  is  satisfied  when  the  phase 
space  is  a  cylinder  whose  axis  coincides  with  the  y-axis. 


In  the  case  of  a  cylindrical  phase  space  there  may  occur  two 


' 
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kinds  of  closed  paths. 

(a)  Closed  paths  of  first  kind.  These  are  the  ordinary  closed 
paths  surrounding  the  equilibrium  states  without  going  around  the 
cylinder  and  they  are  completely  analogous  to  the  closed  paths  in  the 
phase  plane. 

(b)  Closed  paths  of  second  kind.  They  go  around  the  cylinder  with¬ 
out  surrounding  the  equilibrium  states. 

Obviously  both  kinds  of  closed  paths  correspond  to  periodic 

motion. 


7 .  Case  p(t)  -  p  a  constant. 

The  equation  with  which  we  are  dealing  now  is 


x '  =  y 

y'  =  -  cy  -  a  sin(x+p)  +  a  sin  p  . 


The  isocline  (the  locus  of  points  where  the  paths  have  a  given  slope) 
=  0  has  the  equation 


y 


-a  sin(x+p)+q  sin  p 

c 


(4.14) 


and  represents  a  displaced  sinusoidol. 


It  crosses  the  x-axis  only  when 


dy 

a  sin  p  <  a  ,  that  is,  sin  p  <  1.  Also  >  0 

and  and  the  x  axis  and  4^-  <  0  elsewhere  for 

dx 


between  the  sinusoid 


iZ  >  0  if  ~cy  -  a  sin(x+p)  +  a  sin  p  >  Q 
dx  y 


y  >  0 


y 


a  sin  p  -  a  sin(x+p) 

c 


y  <  0 


y  > 


a  sin  p  -  a  sin(x+p) 

c 


The  singular  points  are  at  sin(x+p)  =  sin  p,  y  =  0  and  as  we  saw  before 
are  the  points 


Figure  4.6.  Displaced  Sinusoid. 
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We  consider  first  the  points  B^.  Let  x  =  2mr  +  n  (n  >  0  and  very 
small)  and  develop  sin  (x+p)  in  powers  of  n  •  If  we  keep  only 
terms  of  first  order  in  n,  the  behaviour  of  the  system  in  the 
neighbourhood  of  is  given  by 

dy  _  -cy  -  an  cos  p  +  .  .  . 
dn  y 

Thus  the  coefficient  matrix  is 

0  1 
-  a  cos  p  -c  i 

2  v 

We  now  have  that  A  =  c  -  4a  cos  p,  p  =  -  c,  q=a  cos  p  . 

From  the  table  in  Chapter  I,  we  have,  a  node  if  a  cos  p  >  0 
2  2 

and  c  -  4a  cos  p  >  0,  c  >  4a  cos  p.  The  node  is  stable  for 
%  2 

p  =  -c  <  0.  We  have  a  stable  spiral  point  if  c  <  4a  cos  p.  Note 
when  c  =  0,  we  have  a  centre  as  was  seen  previously  when  we 
examined  the  equation  with  no  damping. 

We  consider  now  the  singular  points  :  x  =  nir  -  2p, 

y  =  0.  Developing  sin(x+p)  in  powers  of  r)  ,  we  get  as  above 

dy  _  -cy  +  an  cos  p  +  . . . 
dn  y 


In  this  case  we  have  a  saddle  point.  Again  when  c  =  0,  we  get  a  saddle 


'a 
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point  as  in  the  case  with  no  damping. 

Andronow  and  Chaikin  [1]  discusses  in  detail  the  existence 
of  periodic  solution  for  such  a  system. 

8 .  The  case  p(t)  a  nonconstant  function  of  t- 

x*  =  y 

y»  =  _  Cy  _  2a  ( cos  sin  y)  . 

We  use  the  same  Lyapunov  function  as  given  by  equation  (4.10)  for 
in  this  case 

V'(x,y,t)  =  -  cy2  -  xYg*y  ‘0S  P(t)  +  “P'  WslnjCt)) 

2  5  4 

+  ax  y  sin  p(t)  -  ay(cos  p(t)o(x  )+sin  p(t)o(x  )). 

2 

The  only  new  term  is  -cy  and  this  does  not  affect  the  sign  of  V' 
for  it  is  always  negative,  so  we  again  have  stability.  Thus 
Theorem  4.11  goes  through  in  this  case. 

9 .  The  case  when  x  =  p(t,e)  is  a  solution 

The  last  case  to  be  discussed  is 


x"  +  ex'  +  a  sin  x  =  ef(t). 


(4.3) 
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Here  we  let  x  =  p(t,e)  be  a  solution.  We  suppose  e  >  0,  small, 
and  that  p(t,e)  admits  the  expansion, 

p(t,e)  =  PQ(t)  +  ep1(t)  +  e2p2(t)  +  ...  =  p  (t)  +  ep^(t)  +  o(e) . 

After  transforming  we  find  that  equation  (4.3)  can  be  written  as  the 
system 


x '  =  y 

,  _  (x+2p(t,e))  .x  //  icn 

y  =  -  cy  -  2a  cos  - S;  ■  ■  sin  —  .  (4.15) 

X  X 

Now  y '  =  -  cy  -  2a [cos  p(t,e)  -  —  sin  p(t,e)  +  o(x) ]  [—  +  o(x) ]  for 

XX  X 

cos^  +  p(t,e))  =  cos  j  cos  p(t,e)  -  sin  -j  sin  p(t,e) 

2 

=  (1  -  \  +  o(x2))cos  p(t,e)  -  (|  -  (|)3+...)sin  p(t,e) 

X 

=  cos  p(t,e)  -  s^-n  p(t,e)  +  o(x)  . 

We  are  considering  here  only  the  first  terms  in  the  expansions  of  sin  x 
and  cos  x 


y’  =  -cy  -  2a[cos(po(t)+ep1(t)  +  o(e) )  -  ^  sin(po(t)+ep1(t)  +  o(e)  +  o(x)  ]  [|-  +o(x)  ] 

X  X 

=  -cy  -  2a[cos  p  (t)-ep.  (t) sin  p  (t)-  sin  p  (t)+  -j  o(e) 

OX  O^-O^ 


' 
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J  P1(t)  cos  P0(t)+  (x))[^+o(x)]  . 


The  last  line  follows  because,  if  €  is  sufficiently  small ^ 


sin(pQ(t)+ep^(t)+o (e))  *  sin  po(t)cos(ep^(t)+o(e))  +  cos  po(t)sin(ep^(t)+o(e)) 

■  sin  p  (t)(l-o(e))  +  ep1(t)cos  p  (t). 
o  1  o 


Similarly 


cos(po(t)+ep1(t)-to  (e))  =  cos  po(t)cos(ep1(t)+o(e))  -  sin  pQ(t)  sin(ep1(t)+o(e) ) 

■  cos  po(t)(l-o(e))  -  ep1(t)sin  pQ(t) 

2  2 

y’  =  -cy  -  2a[~-  cos  pQ(t)-e  ^  P1(t)sin  p  (t)  -  ~  sin  p  (t)-e  p1(t)cos  pQ(t) 

2  o 

+  o (e)+o(x  )  ] 

y'  =  -cy  -  ax(cos  p  (t)-  j  sin  pQ(t))  +  eaxpa (t)(sin  pQ(t)+  ~  cos  pQ(t)) 

-  a  ~  o(e)  -  2ao (x2) . 


The  problem  now  is  to  find  a  suitable  Lyapunov  function  for  this  system. 
Consider 


V(x,y,t) 


1 

2 


y  + 


_a 

2 


cos  PQ(t)  -  eap1(t)sin  p  (t)  ~ 


(4.16) 


i*  t  r.wo  a  l  J enl  srfT 


' 
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Choose  p  (t)  such  that 
*o 


0  <  6  <_  sin  PQ(t)  <_  1  -  6  ’  0  <  6  <  1,  0  <  6'  <  1 


i.e. 


0  <  arc  sin  6  <_  PQ(t)  arc  sin(l-6’) 


0  <  cos  arc  sin(l-6')  £  cos  pQ(t)  £  arc  sin  6  <  1 
0  <  /l  -  (1-6  '  )2  <_  cos  PQ(t)  /l  -  6^  . 


Let 


0  £  p-^t)  <  M,  pj^(t)  >_  0,  p' (t)  _>  0 


Then  V(x,y,t)  _>  jj-  y^  +  (j  /2  6’  -6’2  -  eaM(l-5  '  ))-y-  3  •  We  want  W(x,y) 
to  be  positive  definite  so  it  is  necessary  that 


a/26 ’  -  6’2  -  eaM(l-6 ’ )  >  0 


M  < 


/26  ’  -  6  ’ 2 

£ ( 1—6 ’ ) 


It  is  now  clear  that 

(i)  V(x,y,t)  is  defined  for  all  t  >_  0 

(ii)  V(0,0,t)  =  0  for  t  >_  0 

,  ,  , _ „  2 

(iii)  V(x,y,t)  W(x,y)  =  |  y  +  a  (/2<5 1  -  <5  'z  -e  MCl-S'))-^"  • 


, 
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It  remain  to  show  that  V’(x,y,t)  <  0  . 


V'(x,y,t)  =  yy’  +  axx' 


, .  N  a  2 
cos  PQ(t)  -  -  X 


p^(t)sin  pQ(t) 


2 

"  y  {p1(t)p^(t)cos  pQ(t)+p^(t)sin  p  (t)} 


-  eap^(t)sin  po(t)xx' 


-  ov  xy  cos  PQ(t) 


2 

X 

+  a  ~  y  sin  p  (t)  +  eaxy  p, (t)  sin  p  (t) 
z  o  1  o 


x  x  2 

+  ea-- yp1(t)cos  pQ(t)  -  a  -  y  o(e)  -  2ayo(x  ) 

a  2 

+  axy  cos  p  (t)  -rx  p'(t)sin  p  (t)  -  eaxyp1(t)sin  p  (t) 

O  Z  O  O  X  o 


2 

-  ea  {p1(t)p^(t)cos  PQ(t)+p^(t) sin  p  (t)} 


o  o  ^ 

=  -cy  -  f  X  p^(t)sin  PQ(t)  -eay  {p1(t)p^(t)cos  pQ(t)+p^(t) sin  pQ(t)} 
2  2 

x  x  x 

+  a  —  y  sin  PQ(t)  +  earyyp1( t) cos  pQ(t)  -  a  j  yo(e) 


-  2ayo(x  ) 


We  will  consider  the  last  four  terms  as  higher  order  terms.  The  first 
four  terms  in  the  expression  for  V'  are  negative  so  we  do  not  have  to 


worry  about  them. 


- 
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Since  our  only  concern  is  stability  in  a  small  region  around 

the  origin  say  |x|  <  e^,  |y|  <  the  last  four  terms  are  of  the 

3 

form  c^e  which  can  be  neglected  when  compared  with  terms  of  the  form 
2 

c^e  (the  first  four  terms)  provided  e  is  sufficiently  small.  We 
thus  have  stability  at  the  origin. 

Theorem  4.17  V(x,y,t)  given  by  (4.16)  is  a  Lyapunov  function 
for  the  system  (4.15)  provided 


(i) 

V 

O 

<  sin 

i 

* — i 

V  | 

4J 

'w' 

o 

a 

6’, 

ev 

t — 1 

V 

°o 

V 

O 

(ii) 

p;<o 

1  v 

o 

0  <_  px(t)  < 

M, 

p ' (t)  _>  0  . 

(iii) 

E  >  0 

and 

sufficiently 

small . 

.nigi  to  »r5  Jb  svs  zurll 


(3)^q 

Hi  ns  \  -  l'3 
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APPENDIX 

PERTURBATIONS  ANS  STABILITY 

If  we  consider  two  systems  x'  =  f(x)  ...  (1)  and 

x’  =  f(x,e)  and  if  we  have  some  information  about  the  stability  of  the 
first  system,  it  is  of  interest  to  investigate  the  manner  in  which  the 
stability  of  the  second  system  depends  on  e  . 

When  we  made  the  transformation  on  x'  =  y'  +  p’(t)  and  considered 
the  case  p(t)  =  p  ,  a  constant,  we  were  looking  at  the  system  x’  =  F(x,p) 

. . .  (2)  ,  where  p  could  be  considered  as  a  perturbation  parameter. 

Intuitively,  we  would  expect  that  for  sufficiently  small 
p  ,  the  stability  behaviour  of  the  original  system  should  not  be 
changed . 


This  was  in  fact  borne  out  when  we  considered  the  Duffing 
equation.  Certain  small  value  of  p  had  the  effect  of  shifting  the 
paths  but  not  changing  the  singularities. 

We  should  note  now  that  F(x,0)  =  f(x)  and  F(0,p)  =  0  . 

If  (1)  has  a  solution  cf>(x)  and  (2)  has  a  solution  $(x,p)  then 

$(x,0)  =  4>(x)  since  F  is  in  x  and  p  . 

Let  us  now  examine  x'  =  f(x)  with  the  equilibrium  state 

at  the  origin.  The  perturbed  system  is  x'  =  F(x,p)  .  Let  Q  be 

a  closed  and  bounded  set  containing  the  origin  and  a  subset  of 

Q  .  Let  $  be  a  solution  of  (2)  and  P  be  the  set  of  perturbations 


. 


' 

■ 
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such  that  |p I  <  6  (if  we  were  considering  the  non-autonomous  case 

we  would  want  |p(x,t)|  <  6  for  all  t  _>  0  and  all  x  ).  Suppose 

solutions  which  start  in  Q  remain  in  Q  . 

o 

The  concept  of  stability  is  relative  to  the  sets  Q  ,  Qq 
and  the  number  6  .  Suppose  we  consider  Qq  as  the  original  ball  in 
which  we  have  stability,  we  see  that  6  could  determine  how  close  the 
stability  region  of  the  perturbed  system  would  be  to  Qq  (i.e.  whether 
it  would  be  a  superset  of  Qq  or  be  a  subset  of  it). 

This  question  leads  to  uniform  boundedness  of  solutions  relative 

to  Q  and  the  class  of  perturbations  P  .  The  bound  should  not  only 
o 

exist  but  be  sufficiently  small.  Thus  for  practical  results  a  knowledge 

of  the  size  of  Q  and  6  should  be  known. 

o 

What  this  all  means  is  that  for  very  small  6  ,  we  would 
find,  the  stability  region  is  not  appreciably  changed.  This  has  been 


noticed  in  the  cases  dealt  with  in  this  thesis. 


i 


t 


